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I.  INTRODUCTION 


Several  fast  aerodynamic  prediction  codes  for  missiles  were  written  in 
the  last  decade.  These  codes  were  Intended  for  aiding  missile  designers  In 
obtaining  quick  engineering  estimates  for  the  aerodynamic  coefficients  and  the 
dynamic  stability  of  their  particular  configurations. 

These  fast  codes  are  based  on  1)  basic  and  simplified  theorems,  2)  exper¬ 
imental  data  which  are  algebraically  or  numerically  fitted,  and  3)  some 
empirical  formulae  based  on  observations  and  some  personal  experience.  Metho¬ 
dology  of  these  codes  Is  based  on  missile  component  build-up  with  adjustments 
for  component  interference  (interaction)  effects.  These  codes  were  required 
to  be  fast,  usually  using  less  than  60  CPU  seconds  on  a  typical  mini-computer 
(such  as  a  VAX-11/780)  for  each  flight  condition.  They  originally  were  meant 
to  give  estimates  for  the  basic  aerodynamic  coefficients,  in  particular: 
Cq,  Cn,  Cm,  Cm  ,  CN  over  a  range  of  Mach  numbers  and  angles  of  attack. 

a  Cl 

Now  It  Is  required  that  these  codes  yield  more  accurate  predictions,  to 
provide  all  the  aerodynamic  coefficients,  and  to  cover  a  larger  variety  of 
missile  configurations.  It  is  also  necessary  to  examine  the  application  of 
such  missile  codes  to  gun-launched  projectiles,  both  for  spin-  and  fin- 
stabilized  configurations.  For  this  application  the  L/D  ratio  is  usually 
smaller  than  those  of  missiles. 

At  present,  due  to  more  sophisticated  projectile  and  missile  appli¬ 
cations,  there  Is  a  desire  to  develop  such  codes  to  provide  more  accurate 
predictions,  rather  than  merely  a  rough  tool  to  yield  engineering  estimates. 
To  be  useful  in  that  sense,  the  following  accuracy  guidelines  for  tire  basic 
five  coefficients,  should  be  targeted. 


Cq  within  +  5“ 


CN  ,  CM  within  +  15% 
a  n 


C5  and  (CM  +  within  +  25% 


These  demands  of  accuracy  are  more  relaxed  than  the  accuracy  achieved  in 
actual  firing  tests  In  the  ranges  as  provided  by  Rogers.1  This  relaxation  is 
Intentionally  allowed  because  codes  cover  a  wide  variety  of  body  configur¬ 
ations  «nd  different  speed  regimes  where  different  methods  mey  be  used  and 
extrapolation  of  experimental  data  may  be  allowed.  Rogers1  estimates  the 
accuracy  of  the  free-flight  measurements  to  be  as  follows: 


i 


1 


within  +  1% 


CN  within  +  5% 

a 

CM  within  +  2% 

a 

(C  ♦  C  )  within  +  15% 
q  4 


This  required  accuracy  of  the  codes  has  not  been  achieved  for  the  present 
application  as  will  be  discussed  In  section  IV  for  the  results.  However,  for 
more  traditional  configurations  In  the  low  supersonic  speed  regime  (Mach  num¬ 
ber  1,5  to  2.5),  the  results  are  usually  more  accurate  and  can  fall  within 
these  targeted  accuracy  guide  lines. 

It  Is  the  purpose  of  this  work  to  gauge  the  results  of  the  two  fast  codes 
based  on  the  results  obtained  through  an  application  to  the  hybrid  missile- 
projectile  configuration  of  the  Copperhead.  The  Copperhead  projectile  is  a 
laser  guided,  gun-launched  projectile  with  two  sets  of  spring-out  fins.  The 
geometry  will  be  discussed  in  detail  in  the  next  section. 

The  two  codes  examined  are  the  Naval  Surface  Weapons  Center  Aerodynamic 
Prediction  (NSWCAP)  Code2  and  the  Air  Force  Missile  DATCOM  Code.3  The  former 
code  was  developed  during  the  70's  and  provided  a  good  tool  for  design  config¬ 
uration  studies.  The  latter  code  is  a  more  recent  code  which  is  built  to  make 
use  of  all  the  methodologies  of  the  former  code,  with  modifications  and  im¬ 
provements.  The  code  was  built  to  reflect  updated  theories,  include  more 
recent  and  accurate  experimental  data,  add  more  options  for  practical  missile 
applications  (such  as  non-axi symmetric  bodies,  effects  of  inlets  and  rocket 
motor  thrust),  and  reconstruct  the  code  into  a  more  modular  form. 

Several  studies  were  made  by  different  researchers  with  regard  to  the 
capabilities  of  several  existing  fast  prediction  codes.  Some  of  these  codes 
have  narrow  capabilities  in  terms  of  applicable  configurations,  flight  speed, 
estimating  specific  coefficients  only,  amoiiy  other  restrictions.  Reference  4 
compares  the  capabilities  and  results  obtained  using  MISSILE-2  and  DEMON- 
Series  codes.  Reference  5  list,  and  compares  some  of  the  methods  In  ten 
different  codes,  among  them  are  the  NSWCAP  and  Missile  DATCOM  codes.  Refer¬ 
ence  6  evaluates  the  NSWCAP  and  MISSILE-2  codes  and  refers  to  several  other 
codes.  Reference  7  evaluates  methods  used  for  component  build-up  that  were 
later  used  In  the  Missile  DATCOM  code.  Reference  8  Is  a  description  of  the 
NSWCAP  code,  Its  capabilities,  and  Its  analytical  techniques  as  viewed  by  Its 
authors.  Reference  9  Is  a  description  of  Missile  DATCOM  code  with  regard  to 
Its  different  methods,  as  viewed  by  its  principal  authors. 

It  is  not  the  purpose  of  this  work  to  survey  or  compare  such  variety  of 
codes  but  rather  to  apply  two  particular  codes,  which  are  of  more  general 
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nature  and  which  are  of  Interest  to  the  Army,  to  a  particular  hybrid  pro¬ 
jectile-missile  configuration.  The  objective  Is  to  assess  the  accuracy  of 
these  two  codes  as  applied  to  this  configuration.  A  second  objective  is  to 
Identify  areas  of  needed  development  In  both  codes  for  possible  future 
Improvements. 


II.  GEOMETRY  OF  THE  COPPERHEAD  PROJECTILE 

The  Copperhead  projectile.  Figure  1,  has  a  total  length  of  54  Inches 
(1371.6  mm)  and  a  diameter  of  6.09  inches  (155  mm).  It  has  a  spherical  nose 
cap  and  a  conical  section  of  semi-vertex  angle  of  12.5°  connecting  the  nose 
and  the  body  sections.  There  Is  an  obturator  ring  at  the  end  of  the  body. 
The  base  of  the  projectile  Is  solid  with  no  holes  In  It. 

The  projectile  Is  laser  guided  with  two  sets  of  spring-out  fins.  The 
rear  fins  (tall)  spring  out  shortly  after  the  projectile  leaves  the  gun 
tube.  The  projectile  travels  In  this  configuration,  usually  called  the  launch 
configuration,  unguided  and  with  a  speed  decreasing  from  Mach  number  of  1.8  to 
about  0.95. 

The  front  set  of  fins  (wings)  springs  out  In  the  subsonic  Mach  range  from 
0.95  -  0.80,  and  the  rear  control  fins  (tall)  are  then  activated  to  guide  the 
projectile  to  Its  target.  The  projectile  Is  said  to  be  In  its  maneuvering 
configuration  at  this  Mach  range  with  both  wing  and  tall  fin  sets  deployed. 

The  rear  fin  geometry  Is  shown  In  Figure  2.  The  fin  is  swept  back  20° 
and  Is  tapered  In  thickness  from  the  root  to  the  tip  section.  The  cross- 
section  near  the  root  Is  of  diamond  shape  with  leading  and  trailing  edge 
rounding.  The  fins  are  controlled  through  stems,  with  0.2  Inches  (5.98  mm) 
clearance  between  the  body  and  the  fir  root.  The  pitching  panels,  fins  number 
2  and  4  of  Figure  3,  are  located  ,75  Inches  (19.05  mm)  ahead  of  the  yaw  fins, 
fins  number  1  and  3  . 

The  front  fins  (wings)  are  similar  to  the  tail  fins  except  for  two  dif¬ 
ferences.  First,  the  semi-span  length  is  7.149  Inches  (181.6  mm),  compared  to 
5.974  Inches  (147.2  mm)  for  the  tall  fins.  Second,  there  Is  no  noticeable 
clearance  between  the  fin  root  section  and  the  projectile  body  surface,  since 
these  fins  are  fixed  and  are  not  used  to  guide  the  projectile. 

Both  sets  of  fins  have  slightly  different  shapes  of  slots  In  the  projec¬ 
tile  body  where  they  are  housed  before  deployment.  Both  sets  of  fins  have 
tip  notches  for  the  releasing  mechanisms  to  hold  the  fins  before  they  are 
sprung  out  from  their  housing  locations.  Geometry  of  both  sets  of  fins  Is 
listed  In  Table  (1). 


III.  APPLICATION  OF  THE  TH0  CODES 

Both  codes  were  applied  for  sea  level  conditions  with  a  Reynolds  number 
of  6.18  X  106  per  Mach  number  per  foot.  For  M  ■  1.8,  the  Reynolds  number 
Is  11  X  106  per  foot. 
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Both  codes  were  applied  for  both  launch  and  maneuvering  configurations  in 
the  range  of  Mach  number  0.3  <  M  <  1.8.  Some  modifications  In  the  fin  geo- 
mmetry  had  to  be  made  to  suit  the  Input  capability  of  each  code.  For  example, 
the  fin  swept  tip  chord  had  to  be  made  horizontal  and  the  semi-span  was  ad¬ 
justed  to  account  for  that.  Also,  the  tall  fin  body  gaps  were  not  considered, 
and  the  tall  fins  were  assumed  to  extend  continuously  to  the  root  section. 
Also,  the  details  of  the  obturator  were  ignored  and  the  obturator  was  modeled 
as  if  It  was  a  small  "bump"  on  the  body,  with  a  certain  height  as  Is  usually 
the  case  for  simulating  a  “rotating  band". 

The  zero  lift  case  was  always  computed  In  addition  to  the  small  angle  of 
attack  case  (a  -  2®). 


IV.  RESULTS  AND  COMPARISONS 

Free-fllght  data  are  available  In  Reference  10,  while  wind  tunnel  results 
are  obtained  from  Reference  11. 


1.  LAUNCH  CONFIGURATION 

First,  four  flight  conditions  were  chosen  from  Reference  10  and  both 
codes  were  run  at  Mach  number  and  angle  of  attack  of  (1.77,  2.9°),  (1.47, 
1.8°),  (1.20,  1.1°)  and  (0.81,  0.9°).  The  results  for  Cn,  CK  ,  CN  and  Xc 

o  a  P 

are  considered  reasonable.  The  results  of  (C*  +  CM.)  as  obtained  by  NSWCAP 

q  a 

Is  largely  Inaccurate  especially  for  M  »  0.8  .  For  subsonic  speeds,  the 
NSWCAP  code  does  not  Include  CM>,  therefore  the  value  of  (C^  +  CH.)  Is  not 

a  q  a 

properly  calculated  In  that  speed  regime.  In  fact,  for  the  case  of  (M  »  0.81, 
a  «  0.9')  the  range  result  showed  an  unstable  flight  condition  based  on  pitch 
damping,  while  the  code  predicts  a  stable  condition.  Range  data  are  compared 
to  the  computed  results  In  Table  (2). 

Second,  the  two  codes  were  applied  In  the  Mach  range  of  0.3  to  1.8  and  at 
zero  angle  of  attack.  The  results  for  Cq  Is  shown  in  Figure  4.  Both  codes 

underpredict  the  wind  tunnel  and  range  data.  This  may  be  expected  due  to  lack 
of  consideration  of  the  effects  of  the  fin  slots  of  the  projectile  body.  In 
both  codes.  Also  body-fin  clearance  (gap)  effects  which  should  be  applied  to 
the  tall  fins  are  not  considered  by  either  code.  In  addition,  the  DATCOH  code 
does  not  Include  the  obturator  effect,  which  is  usually  modeled  as  a  rotating 
bend.  The  computed  results  of  both  codes  agree  better  with  the  experimental 
data  In  the  supersonic  regime  (H  >  1.2),  they  worsen  In  the  transonic  regime 
(M  «  0.8  to  1.2)  and  they  deteriorate  further  at  subsonic  speeds  (M  <  .8). 

Reference  12  was  first  to  report  the  effects  of  fin  slots  on  the  normal 
and  axial  forces  of  the  Copperhead.  Wind  Tunnel  tests  were  made  on  a  full- 
scale  projectile  at  both  subsonic  (M  »  0.8)  and  supersonic  (M  -  1.5)  speeds. 

References  13-16  have  also  reported  the  effects  of  body  slots.  Such 

Information  should  be  used  In  the  future  for  modeling  in  both  codes.  Also, 
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Reference  15  suggests  a  modification  to  account  for  the  fin-body  gap  (clear¬ 
ance)  effects. 

Figure  5  shows  the  slope  of  the  normal  force,  Cjj  ,  as  it  varies  with 

a 

the  Mach  number.  The  two  codes  gave  close  values  to  each  other  but  they  both 
considerably  overpredicted  the  range  results  In  the  transonic  regime  between 
Mach  number  0.8  and  1.2.  It  is  surprising  that  the  wind  tunnel  results  are 
also  significantly  higher  than  those  of  the  free  flight  range  tests.  The 
normal  force  predictions  of  the  codes  can  be  improved  if  the  fin  gap  effect 
has  been  accounted  for  and  if  an  average  roll  orientation  angle  is  considered. 

Figure  6  shows  the  slope  of  the  pitching  moment  about  the  C.S. 

Consistent  with  the  overprediction  of  CN  ,  both  codes  overpredict  the  pitching 

(1 

moment  slope.  The  predictions  are  twice  or  three  times  larger  than  free- 
f light  data.  The  OATCOM  code  is  closer  to  the  experimental  data  than  the 
NSWCAP,  due  to  better  prediction  of  the  location  of  the  center  of  pressure. 
The  same  dilemma  of  the  wind  tunnel  data  being  considerably  higher  than  the 
range  data  Is  also  observed. 

Figure  7  shows  the  DATCQM  results  for  the  Xr  location  to  be  more 

LP 

accurate  then  those  of  NSWCAP.  Compared  to  free-fllght  data  of  Reference  10, 
the  OATCOM  results  are  more  accurate,  but  still  overpredict  Xr  by  about  0.4 
calibers.  UP 

Figure  8  shows  the  NSWCAP  predictions  for  the  pitch  damping  coefficient. 

The  OATCOM  code,  on  the  other  hand,  does  not  compute  this  derivative.  The 

trend  shown  agrees  with  the  range  results  only  In  the  supersonic  Mach  range 
down  to  H  «  1.2  .  The  numerical  values  are  about  67*.  larger  than  those 
measured  In  the  free-fllght  range.  It  Is  suggested  that  the  unsteady  pitch 
damping  coefficient,  CM, ,  Is  largely  In  error  possibly  due  to  fin  flutter  or 
a 

to  unsteady  flow  effects  in  and  out  of  the  body  slots  and  arcund  fin-body  gaps 
which  are  not  considered  in  the  code.  Huweve'-,  for  transonic  and  subsonic 
speeds,  the  code  falls  to  predict  the  trend  as  well  as  the  values.  The  lack 
of  Including  C..  for  those  speed  regimes  Is  a  possible  reason  for  such 
failure.  no 


2.  MANEUVERING  CONFIGURATION 

With  the  wing  fins  deployed,  the  projectile  decelerates  from  Mach  0.95 
down  to  Mach  0.3.  Computations  were  made,  however,  for  this  B-W-T  configur¬ 
ation  for  the  Mach  range  of  1.8  to  0.3. 

Figure  9  shows  the  total  drag  coefficient  for  this  configuration  in  com¬ 
parison  to  the  launch  configuration  (B-T).  The  Increase  In  drag  Is  due  to  wing 
fin  drag  less  the  reduction  In  drag  due  to  the  Interference  of  the  wing  fins 
on  the  t«H  fins.  The  DATC0M  code  shows  smaller  Increase  than  that  of  the 
NSWCAP  code,  due  to  vortex  tracking  corrections  included  In  DATC0M,  while  the 
larger  effect  as  computed  by  the  NSWCAP,  Is  due  to  the  lack  of  consideration 
of  wing-tall  Interference  effects.  It  should  be  pointed  out  that  a  recent 
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nonlinear  vortex  tracking  procedure  has  been  developed17  and  proved  to  give 
more  accurate  predictions. 

Figure  (10a)  shows  the  normal  force  slope,  where  the  increase  caused  by 
the  wing  fin  lift  Is  smaller  for  the  DATCOM  code  than  the  increase  predicted 
by  the  NSWCAP  code.  The  cause  for  this  Is  the  reduction  in  lift  of  the  tail 
fin  due  to  the  trailing  vortex  of  the  wing,  as  accounted  for  in  the  DATCOM 
code.  Figure  (10b)  shows  the  change  in  normal  force  slope  as  predicted  by 
DATCOM  Code,  due  to  the  deployment  of  the  wing  fins.  Three  Mach  numbers  .95, 
.9  and  .8  were  chosen  for  the  projectile  speed  at  deployment. 


Figure  (11a)  shows  the  pitching  moment  slope  for  a  range  of  Mach  num¬ 
bers.  For  the  B-W-T  configuration,  the  wing  normal  force  pushes  the  center  o1 
pressure  forward  towards  the  nose,  thus  causing  the  pitching  moment  about  the 
C.G.  to  be  smaller.  Thus  the  projectile  Is  less  stable.  Figure  (lib)  dis¬ 
plays  a  decrease  in  the  dynamic  stability  of  the  projectile  due  to  the  reduc¬ 
tion  in  pitching  moment  slope  from  -25.  to  -5.  The  location  of  the  center  Of 
pressure,  Xc  ,  is  shown  in  Figure  (12a)  to  shift  towards  the  C.G.  and  away 

from  the  projectile  base.  Figure  (12b)  shows  the  sudden  shift  In  the  location 
of  the  Xr  due  to  wing  deployment. 

LP 


The  dynamic  stability  for  pitch  disturbance  remains  almost  unchanged  for 
the  B-W-T  configuration  (compared  to  the  B-T)  in  the  supersonic  regime  as  pre¬ 
dicted  by  the  NSWCAP  code  and  shown  in  Figure  13.  The  DATCOM  code,  on  the 
other  hand,  does  not  compute  this  derivative.  The  trend  shown  agrees  with  the 
range  results  only  in  the  supersonic  Mach  range  down  to  H  «  1.2.  The 
numerical  values  are  about  67%  larger  than  those  measured  in  the  free-fllght 
range.  It  is  proposed  that  the  unsteady  pitch  damping  coefficient,  C^  ,  is 


largely  in  error  possibly  due  to  fin  flutter  or  to  unsteady  flow  effects  in 
and  out  of  the  body  slots  and  around  fin-body  gaps  which  are  not  considered  In 
the  code.  However,  for  transonic  and  subsonic  speeds,  the  code  falls  to 
predict  the  trend  as  well  as  the  values.  The  lack  of  including  C„  for 
those  speed  regimes  seems  to  be  the  reason  for  such  failure.  a 


Figure  14  shows  the  longltudial  stability  chart  for  small  a's  and 
moderate  deflection  angles,  5,  at  Mach  number  0.5.  It  is  shown  that  the 
NSWCAP  code  overpredicts  both  CM  and  C^  for  all  cases,  more  than  the 

DATCOM  code  does.  For  the  same  a,  the  discrepancy  increases  with  increase 
in  6,  Similar  results  are  also  shown  in  Figure  15  for  Mach  number  0.95.  It 
is  noticed  that  the  discrepancy  increased  for  this  transonic  speed  as  was  no¬ 
ticed  earlier  in  Figures  5  and  6  for  CN  and  C^  .  The  DATCOM  code  shows 

a  a 

better  results  than  the  NSWCAP,  especially  for  large  6  due  to  the  inclusion 
of  the  equivalent  angle  of  attack  approach  of  Reference  T/. 

The  roll  damping  coefficient  was  computed  for  both  configurations  only  by 
the  NSWCAP  code  since  the  DATCOM  code  does  not  presently  have  this  capabi¬ 
lity.  The  results  for  B-T  configuration  are  shown  in  Table  3,  where  reason¬ 
able  agreement  with  the  wind  tunnel  results  can  be  observed  especially  when 
excluding  the  transonic  speed  range.  However,  the  results  become  extremely 
large  for  the  B-W-T  configuration,  and  is  attributed  to  lack  c."  consideration 
of  wing-tall  Interference  effects  In  that  code, 
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V.  AREAS  OF  NEEDED  DEVELOPMENT 


In  Table  4,  a  list  Is  given  for  areas  of  needed  development  In  both 
codes.  This  list  was  compiled  through  the  application  to  the  Copperhead  pro¬ 
jectile  case  as  well  as  to  other  cases.  The  order  in  which  they  are  listed 
does  not  reflect  the  order  of  Importance,  because  the  latter  depends  on  the 
objectives  of  each  user  of  the  codes. 


VI.  CONCLUSIONS 

Through  the  application  of  t.he  two  codes-NSWCAP  and  Missile  DATCOM-to.the 
C’  ,'rhead  projectile  geometry,  the  following  conclusions  have  been  drawn. 

1.  The  DATCOM  code  generally  gave  slightly  better  results,  compared  with 
experiment,  than  those  of  the  NSWCAP. 

2.  Both  codes  badly  estimated  the  slopes  of  the  normal  force  and 
pitching  moment  coefficients  due  to  fin  slot  and  fin  gap  effects  which  are  not 
included  In  either  code. 

3.  The  effects  of  the  deflection  angles  of  the  control  surfaces  are  no., 
explicitly  computed  In  either  code.  Both  codes  failed  to  provide  this  Infor¬ 
mation  which  is  essential  to  guided  projectile  configurations. 

4.  The  dynamic  derivatives  of  the  NSWCAP  code  are  not  accurate  for  this 
configuration.  Furthermore,  they  are  not  calculated  In  the  present  version  of 
the  DATCOM  code, 

5.  Both  codes  gave  poor  estimates  for  all  aerodynamic  coefficients  in 
both  the  subsonic  (M  <  .8)  and  transonic  (0.8  <  M  <  1.2)  speed  regimes. 

6.  The  DATCOM  code,  being  developed  more  recently,  Is  written  in  a 
modular  form  allowing  ease  of  modification  and  checking.  The  NSWCAP,  being  a 
pioneer  code,  lacks  this  feature. 

Other  areas  of  needed  development  in  both  codes  were  identified  and  listed  in 
Table  4  for  future  development,  ihese  codes  serve  an  important  function  and 
should  be  developed  to  better  meet  user's  needs. 
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TABLE  1.  Copperhead  Wing  and  Tail  Fin  Geometry. 


Wing  Fin 

Semi-Span  (exposed)  (in) 

Chord  (root  and  tip,  theoretical 
parallel  to  body  (In) 

Area  (single  panel,  one  surface) 

Sweep  Angle  (degrees),  baseline 

Root  chord  thickness  ratio 

Tip  chord  thickness  ratio 

Leading  edge  location  of  root  chord  (in) 

Dimensions 

7.149  (181.58  mm) 

3.051  (77.49  mm) 

20.309  (1.31  x  10“  mm2) 
20 

0.0743 

0.0197 

32.32  (820.93  mm) 

Tail  Fin 

Semi-span  (exposed)  (In) 

5.974  (151.74  mm) 

Chord  (root  and  tip,  theoretical) 

3.051  (77.49  mm) 

parallel  to  body  (in) 

Area  (single  panel  exposed)  (in2) 

16.891  (1.09  x  10“  mm2) 

Sweep  angle  (degrees),  baseline 

20 

Root  chord  thickness  ratio 

0.0743 

Tip  chord  thickness  ratio 

0.0196 

Leading  edge  location  of  root  chord  (in) 

Fins  1,3 

48.640  (1235.47  mm) 

Fins  2,4 

47.992  (1218.99  mm) 

TABLE  2.  Comparison  of  Code  Results  with  Measured  Data 
Launch  Configuration  (B-T) 


Mach  Number, 
Angle  of  Attack 

Predict Ipn 
Method 

CD 

C\ 

(C.G.j 

Rad 

c«„ 

Rad 

XCP 
(Cal  - 
Base) 

rTcrr” 

9 

Sec/Rad 

M  -  1.77, 

a  -  2.9° 

Range 

Test  Results* 

.740 

-0.06 

5.51 

3.69 

-99 

NSWC  Code 

.654 

-2.930 

7.253 

3.30 

-210.12 

DATC0M  Code 

.698 

-1.77 

6.915 

3.454 

t 

M  *  1.47, 

a  -  1.8° 

Range 

Test  Results* 

.760 

-0.88 

5.07 

3.53 

-200 

NSWC  Code 

.671 

-9.606 

8.445 

2.56 

-228.9 

0ATC0M  Code 

.733 

-6.648 

8.073 

2.89 

— -t 
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Continued 


Continued 


M  -  1.20, 

-  1.1* 

Range 

Test  Results* 

.803 

-10.52 

6.96 

2.24 

-132 

NSWC  Code 

.663 

-22.14 

11.33 

1.74 

-248.6 

DATCOM  Code 

.746 

-18.53 

11.35 

2.08 

— — t 

M  »  0.81 

o  -  0.9" 

Range 

Test  Results* 

.398 

-10.56 

8.31 

2.43 

15 

NSUC  Code 

.296 

-22.78 

11.88 

1.78 

-252.8 

DATCOM  Code 

.320 

-18.14 

10.33 

1.95 

- 1 

*  R.  McCoy,  March  1981,  Reference  10. 
t  DATCOM  Code  does  not  compute  this  coefficient. 


TABLE  3.  NSWCAP  Code  Results  for  Roll  Damping  Cofficlent 

CA  [RAD/SEC}-1 
P 


Launch  Configuration  (B-T) 

«  -  0°,  a  »  0° 

Mach  Number 

0.5 

0.8 

0.9 

0.95 

1.2 

1.5 

NSWCAP  Code 

-14.07 

-15.84 

-19.06 

-20.16 

-20.65 

-14.6 

Wind  Tunnel 

-10.50 

-11 

-11.4 

-12 

-16 

-11.1 

Maneuvering  Configuration  (B-W-T) 

6  -  0",  a  «  0" 

Mach  Number 

0.5 

0.8 

0.9 

0.95 

1.2 

1.5 

1.8 

NSWCAP  Code 

-34.97 

-38.74 

-45.0 

-47.0 

-50.92 

-38.5 

-30.49 

Wind  Tunnel 

-20 

-22 

-23.5 

-24.2 

-28 

-25.2 

-23 

TABLE  4.  Capability  Comparison  and  Areas  of  Needed  Development 


NSWCAP  Code 

Missile  DATCOM  Code 

I.  Fins 

la) 

Only  2  or  4  fin  panels  only.  In 

la) 

Only  2  or  4  fin  panels*. 

cruciform  “Plus"  position  only 

arbitrary  roll  angle 

b) 

No  roll  angle  aerodynamics 

lb) 

Arbitrary  roll  orientation 

2) 

Limited  to  two  sets  of  fins 

2) 

Limited  to  two  sets  of  fins** 

3) 

No  body  fin- slot  effects 

3) 

No  body  fin- slot  effects 

4) 

No  fin-body  9ap  effects 

4) 

No  fin-body  9*P  effects 

5) 

No  fin  side- sweep  angle  effects 

5) 

No  fin  side-sweep  angle 
effects 

6) 

No  Interdigitated  wing  and  tail 

6) 

No  interdigitated  wing  and 

fins 

tall  fins 

7) 

No  aft-body  fins 

7) 

No  aft-body  fins 

8) 

No  wing-tail  Interference 

8) 

Includes  a  linear  vortex 
correction  for  down-wash 
effects 

9) 

No  wrap-around  fins 

9) 

No-wrap  around  fins 

10) 

Limited  fir.  cross-section 

10) 

Limited  fin  cross-section 

geometry  options 

geometry  options 

11) 

Gives  erroneous  results  for 

U) 

Gives  much  worse  results 

perfect  delta  fin  (or  close  to 

for  perfect  delta  fin  (or 

perfect  delta  planform 

cl ose-to- perfect  delta 
planform 

12) 

Only  tip  and  root  fin  cross- 

12) 

Multi  fin  cross-section 

sections  be  specified 

geometries  can  be  specified 
(Max.  of  10) 

13) 

Assumes  parallel  line  of  sources 

13) 

Does  not  assume  parallel 

for  fin  geometry 

line  of  sources  for  fin 
geometry 

14) 

Does  not  Include  lifting  surface 

14) 

Includes  the  equlvelant 

non-linearity  at  high  angle  of 

angle  of  attack  for  non- 

attack 

linearity  at  high  a 

II. 

Body  Aerodynamics 

la) 

Computes  base  pressure  drag 

la) 

Computes  base  pressure  drag 

~-l-«  4a-  ..1.1 

or  drag  forces 

b) 

Base  pressure  drag  deteriorates 

b) 

Base  pressure  drag  is  not 

at  large  a  (>10°) 

\Ov«i  pi  edl  cted) 

function  of  « 

♦Arbitrary  number  of  fins  capability  is 

now 

being  added  to  the  newer  version 

of  the  code. 

**A  third  set  of  fins  Is  being  added  In  the  newer  version  of  the  code. 

***P.  esently  beln.,  added  In  the  newer  version  of  the  code. 

Continued 


2) 

No  surface  roughness  or  grooving 

2) 

Includes  surface  roughness. 

effects 

but  no  grooving  effects 

3) 

Includes  rotating  band  contri- 

3) 

Does  not  Include  rotating 

butlon  to  Op 

band  effect  on  Cp 

4) 

Calculates  high  Mach  number  cases 

4) 

Does  not  accept  ar\y  nose 

for  blunt  nose 

bluntness  at  high  super¬ 
sonic  speeds  (M  *  4-5) 

5) 

Yields  fair  blunt-nose  hypersonic 

5) 

Yields  very  poor  blunt-nose 

aerodynamics  (M  >  3) 

hypersonic  aerodynamics 
(M  >  3) 

6) 

Yields  poor  subsonic  and  tran- 

6) 

Yields  fair  subsonic  and 

sonic  blunt  nose  aerodynamics 

transonic  blunt  nose 

(M  <  1.2) 

aerodynamics 

7) 

No  forebody  vortex  shedding 

7) 

No  forebody  vortex  shedding 

effects 

effects 

8) 

No  intermediate  body  vortex 

8) 

No  Intermediate  body  vortex 

shedding  effects 

shedding  effects 

III. 

Vehicle  {Body  and  Fins)  Dynamics 

1) 

Computes  roll  damping  C.  and 

1) 

Does  not  compute  any  dynamic 

pitch  damping  (CM  +  CM  ) 
coefficients  q 

derivatives*** 

2) 

CM  Is  not  computed  for  subsonic 
“a 

or  transonic  speeds  (M  <  1.2) 

{set  to  zero) 

3) 

(CM  +  CM  )  Is  not  adjusted  to  in- 

q  4 

elude  effects  of  deflection  angle 
of  fins  (i.e.  It  remains  constant 
with  6) 

4) 

C.  Is  fairly  computed  for  one  set 

of  fins  only.  However,  It  Is 
largely  In  error  for  wing-tall 

combination,  (no  wing-tail  Inter¬ 
ference  effects) 

IV. 

1) 

Fin  Control 

Only  two  fins  allowed  pure  pitching 

1) 

Independent  fcur-fln 

No  simultaneous  yawing  or  combined 
yawing/pitching 

2) 

deflection  angles 

2) 

No  expressions  or  derivatives  for 

No  expressions  or  deriva- 

control  surface  effectiveness; 

tlves  for  control  surface 

(cN4*  CMfi) 

effectiveness  (CN  ,  ) 

6  6 

Continued 
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3)  (CH  +  }  Is  not  corrected  for  5 

q  & 

(remains  constant  with  variations 
in  6) 

4)  CH  (  and  CM  )  for  ar\y  case  with  fin 

a  a 

deflection  is  calculated  as  C^/aa, 

and  is  void  when  a  ■  0.0 

3)  (CM  +  )  Is  not 

q  & 

calculated 

4)  No  difficulty  in  computing 

Cj^  and  for 

’a  a 

configurations  with  control 
surface  deflection 

V.  General  Features 

1)  Takes  about  40  CPU  seconds  for  a 

1)  Faster  by  a  factor  of  1.5 

single  Mach  number  and  angle  of 

(approximately) 

attack  case  (on  a  VAX-11/780) 

2)  Accepts  a  single  angle  of  attack, 

2)  Accepts  several  Mach  numbers 

and  performs  a  loop  for  up  to  20 

and  performs  a  loop  for  many 

Mach  numbers 

angles  of  attack  {minimum  of 

two)  for  each  Mach  number 

3)  Has  no  difficulty  with  redundant 

3)  Gives  erroneous  results  if 

input  data 

redundant  (but  consistent) 

input  data  is  given 

4)  Uses  input  In  feet  only  (combined 

4)  Can  use  either  In,  ft,  cm  or 

with  some  input  in  calibers) 

meter  units 

DIMENSIONS  IN  INCHES 
AND  (MM) 


(positive) 


M-OATCOM  CODE 
NSWC  CODE 


M-DATCOM  CODE 
NSWC  CODE 


M-OATCOM  CODE 
NSWC  CODE 
RANGE  TESTS 


O  M-DATCOM  CODE 
□  NSWC  CODE 


BT  BWT 


M-DATCOM  CODE 


BT  BWT 


BT  BWT 


NSWC  CODE 
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LIST  OF  SYMBOLS 


C0  »  Total  drag  coefficient,  drag  force/0. 5pJJ2Sref 

Cjj  «  Roll  moment  coefficient,  rolling  moment/  <UsrefLref  -  positive  If 
clockwise  (viewed  from  rear  looking  forward) 

C.  ■  3C,/Sp 

p 

Cm  ■  Pitching  moment  coefficient,  pitching  moment/q-Sre.fLref  (positive 
when  nose  up) 

CM  -  3CM/3a  (1/Rad) 

o 

Cn  *  Yawing  moment  coefficient,  yawing  moment/q-Sre^Lre^  (positive 
when  nose  to  right) 

CN  ■  Normal  force  coefficient,  normal  force/ q„Sref 

CN  -  acN/3a  (l/Rad) 

a 

Cy  *  Side  force  coefficient 

D  «  Body  diameter 

Dref  =  Body  reference  diameter 

L  *  Body  length 

Lref  *  Reference  length,  usually  the  body  diameter 
■  Free  stream  Mach  number 
P  »  Spin  (roll)  rate  (radian/sec) 

c|  »  Pitching  motion  rate  (radian/sec) 

q^  -  Free  stream  <(ynam1c  pressure,  0.5pU„2„ 

sref  ’  Reference  area»  ir0  ref^4 

t  »T1me 

X.  »  Location  of  center  of  pressure,  measured  from  the  C.G.  towards  the 
CP  base  of  the  projectile 

a  »  Angle  of  attack,  positive  when  producing  a  positive  normal  force, 
degrees 

a,  -  Total  angle  of  attack.  Including  side  slip  angle,  degrees 
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LIST  OF  SYMBOLS  (Continued) 


3<x/3t 


Fin  deflection  angle  -  for  fin  1,2, 3,4:  positive  when  producing  a 
negative  (counter  clock-wise  rolling  moment  (DATCOM  notation) 

-  for  fins  2,4:  positive  when  trailing  edge  Is 
down  (NSWCAP  notation) 

Roll  angle  of  the  body  cross-  section 

Fin  orientation  angle,  measured  clock-wise  from  the  vertical 
line  of  the  aj  plane 
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